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ABSTRACT 

Shock heating and particle acceleration processes are some of the most fundamental physical 
phenomena of plasma physics with countless applications in laboratory physics, space physics, 
and astrophysics. This study is moti yated by previous o bservations of non-thermal heating 



of heavy ions in astrophysical shocks (jKorreck et alj 2004). Here, we focus on shocks driven 



by Interplanetary Coronal Mass Ejections (ICMEs) which heat the solar wind and accelerate 
particles. This study focuses specifically on the heating of heavy ions caused by these shocks. 
Previous studies have focused only on the two dynamically dominant species, H + and He 2+ . This 
study utilizes thermal properties measured by the Solar Wind Ion Composition Spectrometer 
(SWICS) aboard the Advanced Composition Explorer (ACE) spacecraft to examine heavy ion 
heating. This instrument provides data for many heavy ions not previously available for detailed 
study, such as Oxygen (0 6+ , 7+ ), Carbon (C 5+ , C 6+ ), and Iron (Fe 10+ ). The ion heating is 
found to depend critically on the upstream plasma (3, mass per charge rati o of the ion, M/Q, 
and shock magnetic angle, 9sn- Similar to past studies ( Schwartz et al.l ["l988). there i s no stro ni 



depen dence of ion heating on Mach number. The heating mechanism described in iLee fc Wi 
(2000) is examined to explain the observed heating trends in the heavy ion thermal data. 



Subject headings: Sun: coronal mass ejections (CMEs) — acceleration of particles 



1. Introduction 

Shocks constitute a non-linear transition be- 
tween two dynamically different plasma states. At 
this transition, a portion of the bulk flow energy 
of the plasma is transferred into thermal energy of 
its components. A very small fraction (1-10%) of 
the interacting plasma becomes highly energized 
and is injected into particle acceleration processes 
boosting the particles to very high energy (E> 
GeV). There is, however, a lack of understanding 
of the dynamic effects that lead to this heating. 
The heating of an ion species appears to be depen- 
dent on the ion properties as well as on the bulk 
properties o f the shock and the initial state of 
the p la sma ( Shimada fc Hoshinol 120031: ILee fc Wul 



by the interplanetary manifestations of Coronal 
Mass Ejections (CMEs) - Interplanetary CMEs 
(ICMEs). 

ILee fc Wul |2000) discuss three parameters key 
to understanding heating in collisionless shocks: 
Alfvenic Mach number, Ma, magnetic shock an- 
gle, Obu, and upstream plasma /3, the ratio of ther- 
mal to magnetic energy. The dependence of heat- 
ing on these parameters has been examined be- 
fore in solar w ind shocks both observationally and 
theoretically (lOgilvie et all 19801 Zertsaloy et al 



19761: IBerdichevskv et al]ll997l IZhao et al.lll99lh . 



2000; iKorreck et all |2004 IBerdichevskv et al 



19971 ). It is the purpose of this paper to inves- 
tigate the heating of solar wind ions using accessi- 
ble thermal data of heavy ions from shocks driven 
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Work by IBerdichevskv et al.l (|l997l ) examined the 
heating of ions in the solar wind. The heating 
of heavy ions such as 6+ was found to be more 
than mass proportional when compared to the 
heating of protons. This appears to contradict 
heating seen in other shocks such as those around 
supernova remnant s (TKorreck et"aL 2004 ). The 
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ion and proton thermal velocity were similar, i.e. 
energy is distributed mass proportionally. 

The effect of Alfvenic Mach number of the shock 
on heating has been studied observationally in the 
heliosphere. The Alfvenic Mach number is used 
instead of the acoustic Mach number due to the 



magne tic nature of the solar wind. ISchwartz et al 



( 19881 ) used data from crossings of the Earth's bow 
shock to study the heating of ions and electrons. 
Only a weak correlation between increasing Mach 
number and decreased heating was found. 

The plasma (3 is a d efining characteristic o f 
the plasma. For example, Lipatov fc Zank ( 19991 ) 
found for plasmas with (3 > 0.1 ions were unable 
to "surf" waves as an acceleration process, making 
plasma /3 a determining factor in particle acceler- 
ation models. Note that in t he heliosphere j3 is 
normally < 1.0 (|Galvin Ill997h . 

The magnetic angle from the shock normal, 
9bu, has also been examined with respect to the 
acceleration of ions. Quasi-parallel shocks have 
been shown to accelerate ions from thermal ener- 
gies to well above what w ould be expec t ed from 
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the bulk energy transfer (lEllisonl 119811: ISch olcr 
1990[ iGiacalone et alj Il992t iKucharek &: Scholer 
19951 ). Quasi-perpendicular shocks have been 
shown to be inefficient for accelerating ions from 
thermal ion populations, however these shocks are 
the most efficient accelerators if there is a seed 
population available. The acceleration efficiency, 
once the seed popul ation is available, i s inversely 
proportional to 9bu ( Ellison et al.lll995l ): the more 
perpendicular the shock, the more efficient the ac- 
celeration. Hence, in order to understand cosmic 
ray acceleration, the heating of a suprathermal 
seed particle population is of great importance. 

In addition to the three parameters mentioned 
above, the mass-to-charge ratio also plays a cru- 
cial role in heating of heavy ions. As the mass-to- 
charge ratio incre ases the heating of ions increases 
( Lee &; Wu 2000h . The mass-to-charge ratio dic- 
tates which types of MHD wave interactions are 
available to the particles at the shock front. 

The novel data on ICMEs from the Ad- 
vanced Compositi on Explorer (ACE) spacecraft 
( Stone et alj Il998l ) afford us the opportunity to 
study these shocks and the heating processes for 
individual ion species in greater detail. While the 
previous work included only protons, helium, and 



oxygen (0 6+ ) in the analysis, this study is ex- 
tended to include 7+ , C 5+ , C 6+ , and Fe 10+ ions, 
extending the mass-to-charge ratio range to 1 thru 
5.6. 

This paper examines the heating of heavy ions 
with respect to all critical parameters mentioned 
above that dominate shock heating - the magnetic 
shock geometry - Alfvenic Mach number - 

M^4, plasma j3, and M/Q, the mass per charge 
ratio. We utilize thermal properties measured 
by the Solar Wind Ion Compositio n Spectrome- 
ter (SWICS) (|Gloeckler et al.lll998l ) on the ACE 
spacecraft. We characterize a statistically signif- 
icant sample of 20 shocks for which heavy ion 
heating is observed with high statistical accuracy 
and with well-defined shock properties. Magnetic 
field d ata from the Magn etometer (MAG) instru- 
ment ( Smith et al . 1998) and proton data from 
the Solar Wind Electron Proton Alpha Monitor 
(SWEPAM) instrument (jMcComas et all Il998l ) 
were also used. 

The plasma data available from the ACE space- 
craft and a brief discussion of errors are described 
in Section 2. The selection of ICME shocks is dis- 
cussed in Section 3. In Section 4, the heating of 
ions is examined with respect to each parameter. 
A method of heating is described in Section 5 and 
conclusions are summarized in Section 6. 

2. Observations 



The ACE spacecraft, launched in 1997, orbits 
about the first Lagrangian, LI point, 235 Earth 
radii upstream along the Sun-Earth line. Three 
of the nine instruments aboard were used for this 
study. The SWICS, MAG, and SWEPAM instru- 
ments provide data on ion composition, velocity, 
density, and the ambient magnetic field. For this 
study, the following physical quantities were ana- 
lyzed: proton temperature, proton thermal veloc- 
ity, proton number density, thermal velocities of 
He 2 +, 6 +, 7+ , C 5+ , C 6+ , and Fe 10+ , and mag- 
netic field. 

The heavy ion data (Z>1) was obtained from 
the SWICS instrument. This instrument is com- 
posed of an electrostatic analyzer, which measures 
an ion's energy per charge, and a time-of-flight- 
energy telescope, which measures an ion's velocity 
and total energy. The combination of these two 
sensors enables unique identification of a particle's 



2 



mass, charge, and energy (jGloeckler et al. I ll998h . 
For this study, we use SWICS data with 1-hour 
time resolution; this time resolution provides suf- 
ficient statistics for the minor ions studied. For 
each one hour time bin, five 12-minute accumu- 
lation cycles were summed. Care has been taken 
to exclude the accumulation cycle that contains 
the moment of the shock passage, because it mixes 
upstream and downstream distributions. The pro- 
ton thermal velocity, the proton number density, 
and the magnetic field data were taken from the 
combined SWEPAM/MAG data set available on 
the web at: 

|http://www.srl.caltech.edu/ACE/ASC/| 64-second 
averages of data were used. The data from each 
physical quantity were averaged for one hour be- 
fore or one hour after the shock. This allows for 
a relatively local measurement of the shock, but 
statistically significant heavy ion data. There are 
important statistical limits to these data, which 
are discussed in section 12.11 

2.1. Errors 

For the proposed heating analysis, it is neces- 
sary to quantify the thermal width, or thermal 
speed, of the heavy ion distributions. However, 
the thermal speed is only defined if distributions 
are Gaussian. We have therefore implemented a 
test of the fit of a Gaussian to the data using the 
following technique. For each heavy ion measure- 
ment, the thermal speed was calculated using two 
techniques: by calculating the second moment of 
the observed distribution, v t h, and by determining 
the thermal speed by least-squares fit of a Gaus- 
sian to the distribution, Vthju- The deviations 
between these two quantities are a measure of the 
deviations of the distributions from a Gaussian 
shape. 

The deviations were relatively small for most 
ions with an average value of <10%, but we do 
find a few examples for which energetic tails on 
the distribution function dominate, and the error 
on v t h was then increased, to a more conservative 
error estimate of 20% in order to account for this. 

In order to combine plasma data with our 
composition-resolved dynamics data, we averaged 
the necessary solar wind parameters for 1 hour 
ahead of and behind the shock front using data 
taken at intervals of 64 seconds, to overlap with 
the SWICS data-interval. The statistical error was 



then calculated as the standard deviation of each 
averaged quantity. 

3. Shock Selection 



Sho cks were selected from the lCane fc Richardson 
(|2003l ) list of ICMEs and the ACE Shock List 
maintained online 0. The ACE Shock List details 
the time of the shock passage, the angle between 
upstream magnetic field vector and the shock nor- 
mal, #Bn, and the upstream Mach number, 

The first criterion for selecting a shock for this 
study was to have all crucial data-sets available 
for two hours before and two hours after the shock 
passage. The next criterion for selection was the 
characteristics of the kinetic temperature of the 
solar wind before the shock passage. Since the 
kinetic temperature of the solar wind provides 
a measure of the shock and the plasma charac- 
teristics upstream and downstream of the shock, 
an increase in temperature before the time of 
the shock indic ates pre-heating, ion feedback or 
a rev erse shock (jPaschmann et al.lll981t iGloeckler 
1999). The presence of these features would bias 
the study since ions are heated differently, hence 
these instances need to be excluded. If the tem- 
perature increased in the time period between 60 
to 30 minutes prior to the shock passage by more 
than 30% of the mean value of the temperature 
calculated up to 30 minutes before the shock, the 
shock was considered to exhibit pre-heating and 
was excluded from the analysis. 

Representative solar wind parameters for an 
ICME shock are plotted versus time in Figure [T] 
Panel A shows the solar wind proton velocity, v p , 
versus time as the solid line and the velocity of 
each heavy ion is included as a different symbol, 
as indicated in the plot legend. Panel B shows 
the proton number density, n, in the solar wind. 
Panel C is a plot of the thermal velocity, v t h, of 
protons with the over-plotted symbols represent- 
ing the thermal velocity of individual ions. Panel 
D plots the solar wind proton temperature, T p , 
versus time. Panel E contains the magnitude of 
the magnetic field, B, versus time. Panel F is a 
plot of the magnetic latitude, S, and longitude, A, 
versus time. 



http: //www.bartol .udel . edu/$\sim$chuck/ ace/ACElists/obs_list .html 
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In order to examine the dependence of heat- 
ing on shock orientation with respect to the back- 
ground magnetic field, the shock list was bro- 
ken up into quasi-parallel and quasi-perpendicular 
shocks. Quasi-parallel shocks were defined to have 
&Bn between and 20 degrees and shocks with 9bu 
between 80 and 90 degrees were considered quasi- 
perpendicular. There were 16 quasi-perpendicular 
shocks available for study and 4 quasi-parallel 
shocks. In subsequent sections, each parameter 
will be explored for a perpendicular and parallel 
shock separately. 

4. Heavy Ion Heating 

There are two methods of heating ions that are 
widely invoked in the study of collisionlcss shocks. 
First, the transfer of the bulk fluid kinetic energy 
to the thermal energy of the particles leading to 
mass proportional heating or all ions having the 
same velocity. The second method is wave parti- 
cle interactions that would result in a M/Q depen- 
dence in heating. 

Another point of interest is the distribution of 
heating among the diff erent ion species. In stud - 



ies in the heliosphere (jBerdichevsky et al.l 11997). 



it was found that ions are heated more than mass 
proportionally than ions. However, in studies of 
the collisionless shocks in supernova remnants, 
the ions are hea ted less than mass p ro portionally 



to the protons (jKorreck et all 12004 ) . iLee fc Wul 



(2000) found that the heavy ions are heated by di- 
rect transmission through the shock as reflection 
of heavy ions by the shock potential is prohibitive 
due to their large mass. 

In order to determine which mechanism or 
mechanisms are responsible for the heating, we 
first quantify the heating of each ion species, in- 
cluding protons, and then examine the heating 
dependence on M/Q, M^, and (3 to try to eluci- 
date a candidate mechanism for heating. 

In this study the heating, H, of an ion at a 
shock is defined as the ratio of the square of the 
upstream and downstream thermal velocities. 



H 



v th d 3kT S;d /m a 



v. 



th v 



3&^s,u/ Wis 



^ (1) 



where 

T s = Temperature of the species 
m s = Mass of the species 



v t h d — Average thermal velocity of the species one 
hour downstream 

vth u = Average thermal velocity of the species one 

hour upstream 

k= Boltzmann constant 

However, it has been recognized that the so- 
lar wind is not initially in thermal equilibrium: 
heavy ions hav e a tendency to exhibit high kinetic 
temperatures. dHefti et al. 19981: Zurbuchen et al 



bOOOUvon Steiger & Zurbuchenll2003l) . In order to 
accurately describe the ion heating, the preshock 
temperature ratio of ion temperature to proton 
temperature needs to be considered. The ratio 
of preshock ion to proton temperature shows only 
one of the 96 data points to be significantly less 
than 1. The data in general supports an initial 
condition where ions have equal if not higher ther- 
mal speeds than the protons. Therefore, relative 
heating ratios based on initial thermal speeds will 
give the clearest quantification of heating through 
the shock passage. 

Examining the current data set with respect 
to the temperature ratio of the downstream ions 
to protons, we have found that the heating is 
on average less than mass proportional for the 
ions except Fe +1 °, although 27 out of the 96 data 
points for perpendicular shocks did show greater 
than mass propor tional heating. This contr adicts 
earlier studies by iBerdichevskv et al. (1997) that 
found all of the ion heating to be greater than mass 
proportional. 

The heating of the ions are plotted versus the 
initial temperature ratio of the ion to the proton 
in Figures and ^Bp. The values of initial ratio 
are all greater than or equal to one except for one 
data point. The heating, if greater than one, in- 
dicates that the species was heated through the 
shock. However, only 69% of the ions were heated 
in the perpendicular shocks and 70% in the paral- 
lel shocks. The others actually showed a decrease 
in temperature across the shock front. Protons 
all gained heat or maintained their temperature 
across the shock front. 

In Figures [3^ and [3)3, the ratio of ion and pro- 
ton heating, Hi/H pi is plotted versus the initial ra- 
tio of thermal speeds of the ion and proton where 
Hp is the heating ratio for protons. Only 43% of 
the ions in perpendicular shocks were heated more 
than the protons and 50% of those ions in parallel 
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shocks. The protons are therefore gaining more 
heat than the ions in these shocks. 

To characterize the heating mechanism consis- 
tent with the current data set, the dependence 
of the heating on Alfvenic Mach number, Ma, 
plasma j3, and mass to charge ratio, M/Q, for both 
quasi-parallel and quasi-perpendicular shocks will 
now be examined. 

4.1. Dependence of Heating on Alfvenic 
Mach Number 

Shocks can heat ions by converting bulk kinetic 
energy into thermal kinetic energy of the particles. 
The Mach number gauges the speed of the propa- 
gation of bulk fluid flow relative to the surround- 
ing medium. The Mach number used in this study 
is the Alfvenic Mach number, which takes into 
account interaction with the surr ounding magne- 
tized plasma. Previous studies, ([Schwartz et al.l 
1988) showed the collisionless heating to be weakly 
dependent on Mach number. 

The method of energy transfer can be deter- 
mined by the Critical Mach number. Below an 



Alfve nic Mach number of approximately 2 ([Kennel 
1987), also called a sub-critical shock, the effective 
dissipation of energy by the shock is local, i.e. the 
heating is spatially local to the shock front via 
conduction. Above this critical Mach number the 
dissipation is based on multiple streams of ions 
which have longer characteristic scales. A conve- 
nient measure of a supercritical shock is the ratio 
of the Alfvenic to Magnetosonic Mach numbers. 
If this ratio is > 1 then the shock is supercrit- 
ical by this definition. We have calculated that 
all the shocks in this study are supercritical. The 
data therefore implies a heating method that in- 
volves wave particle interaction. Protons heating 
has been shown to result from both direct trans- 
mission an d relflected prot ons in these supercriti- 
cal shocks ( Lee et al1ll987h . 

Figure [4^, shows the Alfvenic Mach number 
versus the heating for quasi-perpendicular shocks. 
The heating, H, is described by Equation [T] The 
line shows a linear fit to the data with an in- 
creasing trend. The parallel fit seems kinked due 
to the logrithmic scale for the data. The quasi- 
parallel shocks in Figure HJd also shows an increas- 
ing trend with increasing Mach number however 
with a smaller data set and less statistics it seems 



to be a weaker effect. This contradicts the find- 
ings of lSchwartz et al. ( 1988f ). however the fits are 
marginal (x„ ~ 5 for the parallel fit and 0.6 for 
the perpendicular fit). 

4.2. Plasma (3 Effect on Heating 

Plasma (3 is plotted versus the heating for each 
ion species for perpendicular shocks in Figure [S^,. 
The plot shows that with increasing (3 the heat- 
ing of the ions decreases. The heating is therefore 
more effective in magnetically dominated regimes. 
The lines plotted are the fit to an exponential fit 
to the data. 

For the parallel shocks, Figure [5Jd, there is 
an observed decrease in heating with increasing 
(3, however, only a few shocks are available for 
study and hence the trend is more uncertain but 
decreases more sharply than the perpendicular 
shocks. 

4.3. Heating and Mass to Charge Ratio 

The mass to charge ratio of an ion determines 
the type of wave and particle interaction that the 
ion can experience. In Figures [6^ and[6)D, the av- 
erage heating per ion is plotted versus the mass to 
charge ratio for perpendicular and parallel shocks 
respectively. The line is a least squared fit to the 
mean values with the standard deviation of the 
mean values as error bars. For the perpendicu- 
lar shock the fit has a reduced chi-squared value 
of x^=0.25, for the expression H=l. 15+01. 6M/Q. 
The fit for the parallel shocks, xt ~0.48, the line 
is present to emphasize the trend towards an in- 
crease in heating with increasing M/Q ratio. 

In the perpendicular shocks both the M/Q=2 
ions are heated approximately the same, whereas 
in the parallel shocks they are heated differently, 
making the mass important in the parallel shock 
given that the larger value is the C 6+ ion, where 
the M/Q ratio is important for the perpendicular 
shock. 

These shock exhibit an increase in heating 
based on a larger M/Q ratio regardless of the ex- 
act 0Bn, however, the effe ct is more pronou nced in 
the perpendicular shock. iLee fc Wul |2000) found 
that due to transmission of ions across the shock, 
the gyro speed, v g , is approximately equal to the 
thermal speed. The gyro speed increased as the 
M/Q ratio increased which agrees with this data 
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set. 



5. Discussion of Heating Mechanisms 

A shock heating mechanism is needed that ad- 
dresses all observational constraints found in this 
study. Any proposed mechanism must meet the 
following criteria: 

1. Heating increases with decreasing j3 

2. Perpendicular shocks heat more effectively 
than parallel shocks 

3. The ions that initially have thermal veloci- 
ties nearly equal to protons are preferentially 
heated 

4. Heating increases with increasing mass-to- 
charge ratio. 

5. Increasing Mach Number increases the heat- 
ing 



Lee fc Wu| (|2oonh modeled heating in the so- 



lar corona. This heating mechanism and the pre- 
dictions for the downstream thermal speed will 
be used to evaluate this mechanism for the CME 
shock heating. The protons are the main species 
and will form the shock while the heavy ions 
will act as test particles through the shock front. 
The simplifying assumption will be made that all 
changes in the ions velocity will occur parallel to 
the shock normal. The proton speed is also as- 
sumed to be the speed of the solar wind or the up- 
stream bulk speed. As protons pass through the 
shock they are slowed according to the Rankine- 
Hugoniot conditions in order to conserve energy. 
Their bulk kinetic energy is changed into poten- 
tial energy forming an effective electrostatic po- 
tential. Therefore, the potential is proportional to 
the change in energy of the proton. 



1 2 1 2 



P u p,d 



(2) 



This potential jump is seen by the heavy ions as 
they approach the shock front. But the heavy ions 
conserve energy, increasing their velocity across 
the shock according to their M/Q ratio. 

AK.E. Z = q t <j) = ^mi(vl u - (3) 



where is the charge of the ion, raj is the mass 
of the ion in units of proton mass, Vi. u is the veloc- 
ity of the ion upstream and v^d is the velocity of 
the ion downstream of the shock. Plugging in the 
potential formed by the protons © to the kinetic 
energy of the heavy ions ([3]), we can find an ex- 
pression for the downstream ion velocity in terms 
of the M/Q (a) and observed parameters. 



(4) 



where a is the mass to charge ratio, v p is the 
solar wind proton velocity, and c is the ratio of 
downstream proton velocity to upstream proton 
velocity in the shock frame. 

The velocity computed in Equation |4] is the ki- 
netic velocity of the particle at the shock front. 
However, we are trying to examine the thermal 
speed of the heavy ions. As the distribution of 
the ions moves away from Maxwellian, the ther- 
mal speed becom es a measure of kinetic energy. 
Lee fc Wul (l2000h derived that for a perpendicular 
shock the gyro velocity was approximately equal 
to the thermal veloci ty. Using the de finition of 
gyro velocity from the lLee fc Wul ((2000) paper, 



Vi,d,th = Vgyro = V PjU \^ (1 - ^) - (5) 

where a is defined as eA(j)/eo - the ratio of elec- 
tric potential at the shock front to the kinetic en- 



ergy of protons upstream 



is the upstream 



solar wind proton velocity, B t l is the upstream 
magnitude of the tangential magnetic field, and 
B t ,2 is the downstream magnitude of the tangen- 
tial magnetic field. 

By inspection, the downstream thermal veloc- 
ity does increase with increasing M/Q ratio. Us- 
ing the definition of a Mach number the relation- 
ship between the upstream proton velocity and the 
Mach numer is found to be v |)iU =Mav j 4. Insert- 
ing this relation into the square of Equation [5] in 
order to compare with the ion heating data, the 
heating of the ions is proportional to the square 
of the Mach number. However, in the actual data 
sets, there is an increase of heating with respect to 
increasing Mach number however not a significant 
heating as suggested by this equation. 

The plasma (3 is inversely correlated with the 
square of the Alfven velocity. With respect to 



G 



plasma (3 the data match the predicted trend of 
Equation [5] as the heating decreased with increas- 
ing /3. 

As the Mach number increases we expect to 
have greater shock compression of the magnetic 
field which would lead to a larger B ti 2- The Mach 
number is not as strongly correlated as the mag- 
netic characteristics, such as Osn or f3, with the 
heating of the heavy ions, however, the bulk mo- 
tion provided by shocks at higher Mach numbers 
is necessary for the compression of downstream 
magnetic fields that influence the gyro velocity. 

6. Conclusions 

This study examined different ion species in the 
solar wind to investigate the heating that occurs at 
a collisionless shock front ahead of ICMEs. Based 
on the magnetic field angle to the shock normal, 
0Bm each shock was analyzed for its dependence 
on Mach number, M^, plasma (3, and the behavior 
of ions based on their mass to charge ratio, M/Q. 
An increasing [3 was found to reduce the ion heat- 
ing at the shock, highlighting the importance of 
the magnetic field to the heating process. As the 
Mach number increased the heating also increased. 
As the mass to charge ratio increases, the heating 
of the ion increases. Quasi-perpendicular shocks 
were shown to heat the ions more efficiently than 
the quasi-parallel shocks. 

A heating metho d based on a pos t-shock po- 
tential described by iLee fc Wul (|2000h was exam- 
ined for the heating trends observed in this study. 
The data and the predicted values of heating were 
not well correlated for the parallel shock but did 
match trends found for the M/Q, Ma and plasma 
(3 trends for the perpendicular shocks. The mech- 
anism reveals an increased heating with increas- 
ing Mach number. However, the magnitude of the 
increase in heating with Mach number is not of 
the same order of magnitude that is found in the 
data and the trend however wea k is inverse to that 
found bv lSchwartz et al I (|l988l) . 

There is currently similar ion heating data 
available for collisionless shocks in supernova rem- 
nants showing non-preferential heating to ions. 
Using the data compiled here we can compare the 
heating ratios of ions to protons downstream and 
use that ratio to fin d a plasma (3 for the su pernova 
remnant (SNR). In Korreck et al. ( 2004 ). it was 



found that the ratio of oxygen to proton temper- 
ature post-shock was 8.3. Using the fit of heating 
versus plasma (3 for the CME data, we obtain a 
(3 for SN1006 of 1.7. Knowing the temperature, 
T=1.5 xl0 9 K, and the density given, 0.25 cm" 3 , 
we can also find a magnetic field of 1/iGauss. This 
is very close to the assumed galactic magnetic 
field of 3/iGauss. These types of diagnostics could 
lead to a better understanding of the magnetic 
field in SNRs and their role in ion heating and 
acceleration. 

This begs the question of the ubiquitous nature 
of the shock physics: what is the dominant factor 
to determine effect heating at a collisionless shock 
front? The supernova has a Mach number 10 times 
that of the CME shocks however, as seen in this 
data set, the Mach number is not the only fac- 
tor in determining heating. Density and magnetic 
energy seem to be of greater importance. 
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Fig. 1. — Plot of ACE magnetic and temperature 
data versus time, in fraction of a day, for a parallel 
shock. Panel A plots the solar wind proton veloc- 
ity, v p as the solid line and the velocity of each 
heavy ion is included as a different symbol. Panel 
B is plot of the proton number density, n. Panel 
C is a plot of the thermal velocity, v t h , of protons 
with the symbols representing the thermal velocity 
of individual heavy ions. Panel D plots the proton 
temperature, T, versus time. Panel E contains the 9 
magnitude of the magnetic field, B. Panel F is a 
plot of the magnetic latitude, 5, and longitude, A. 
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Fig. 2. — Shock ion heating versus upstream 
thermal temperature ratio, (a) for perpendicular 
shocks and (b) for parallel shocks. The upstream 
ratio of ion thermal velocity to proton thermal ve- 
locity, Tj/Tp, is the x-axis. Ion heating, Hj, is the 
y-axis. 
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(A) Perpendicular Shocks 
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Fig. 3. — Shock heating ratios versus upstream 
thermal temperature ratio, (a) for perpendicular 
shocks and (b) for parallel shocks. The upstream 
ratio of ion temperature to proton temperature, 
Tj/Tp, is the x-axis. The ratio of temperature 
increase between the ion, Hi is the y-axis. 
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Fig. 4. — Plot of Alfvenic Mach number versus 
heating, H, for all the heavy ions. The line indi- 
cates a least squared fit to the data which appears 
kinked due to the logrithmic scale. 
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(A) Perpendicular Shocks 
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(B) Parallel Shocks 
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Fig. 5. — Plot of (3 versus heating, H, for all the 
heavy ions in a shock. The heating is the ratio of 
the square of the downstream ion thermal velocity 
to the square of the upstream ion thermal veloc- 
ity. (3 is the ratio of thermal to magnetic energy 
densities. The line is the exponential fit to the 
data. 
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Fig. 6. — Plot of M/Q versus heating, H, for all 
the ions in a shock. The heating is the ratio of the 
square of the downstream ion thermal velocity to 
the square of the upstream ion thermal velocity. 
The heating for each ion was averaged in order to 
determine a trend in the data. The line is the least 
square fit to the data. 
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